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ABSTRACT 

An  experimental  study  was  carried  out  to  investigate 
the  actual  influence  of  the  boundary  layer  on  heat  transfer 
losses  in  typical  expanders  utilized  in  low  temperature 
systems.   More  specifically,  to  investigate  the  effect  of 
the  boundary  layer  on  the  time  dependence  of  the  rate  of 
heat  transfer  through  the  cylinder  walls.   The  experimental 
expander  utilized  was  a  closed  cylinder  apparatus  with  a 
fixed  mass  charge  of  helium  gas.   Cyclic  pressure  and  volume 
data  points  were  generated  during  operation  of  the  expander. 
The  resulting  experimental  pressure  vs.  volume  data  trace 
was  stored  digitally  and  presented  graphically  using  a 
Nicolet  storage  oscilloscope.   The  CMS  computer  system  was 
used  to  further  reduce  and  interpret  the  data. 

The  results  concluded  that  the  temperature  of  the 
boundary  layer  gas  had  the  controlling  effect  over  the 
cylinder  bulk  gas  temperature  in  determining  the  actual  rate 
of  heat  transfer  through  the  cylinder  walls.   The  apparent 
zero  rate  of  heat  transfer  point  occurred  when  the  bulk  gas 
temperature  equaled  the  cylinder  wall  temperature.   Because 
of  the  strong  effect  that  the  boundary  layer  had  on  the 
heat  transfer  in  the  cylinder,  the  actual  zero  rate  of  heat 
transfer  point  occurred  when  the  boundary  layer  gas  temper- 
ature equaled  the  cylinder  wall  temperature.   The  boundary 
layer  gas  temperature  was  determined  responsible  for  up  to  a 
6k   degree  phase  difference  in  the  times  of  the  actual  and 
apparent  zero  rate  of  heat  transfer  points  in  the  compression 
and  expansion  strokes  in  the  cylinder.   This  phase  difference 
was  evaluated  to  be  leading  in  both  cases.   Therefore  the 
actual  time  of  zero  rate  of  heat  transfer  preceded  the  ap- 
parent time  of  zero  rate  of  heat  transfer  by  a  significant 
margin  in  each  piston  stroke. 
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DEFINITION  OF  SYMBOLS  USFD 

Symbol         Definition  Units 

p 
A  heat  transfer  area  mw 

C  specific  heat  at  constant       3tu/lb  -  R 

"D  m 

pressure 

CT  specific  heat  at  constant       Btu/lb  -  R 

volume 

h  surface  coefficient  of  heat 

transfer 

m  system  mass 

P  cylinder  pressure 

P  reference  cylinder  pressure 

P  compression  stroke  oressure 

C0C1?  at  reference  volume 

P  expansion  stroke  pressure 

*  ^  at  reference  volume 

P  mean  pressure  over  the  com- 

mep  pression  stroke  (eqn  21) 

PSTR  non-dimensionalized  pressure  — 

(eqn  22) 

q  rate  of  heat   transfer  in-lbf/sec 

Q  quantity  of  heat   transferred  in-lb~ 

R  specific   gas   constant  ft-lb~/lb  -°R 

RPM  piston  revolutions  per  minute  rev/min 

S  entropy  of   cylinder  bulk  gas  in-lbf/°R 

SSTR  non-dimensionalized   entropy  — 

(eqn   2k) 

t  time  sec 


Btu/hr-f 

t2- 

•°H 

m 

lb  ./in2 
i 

lbf/in2 

lbf/in2 

lb  J±nZ 

i 

lbf/in2 

Symbol         Definition  Units 

T,  T,  ,,        cylinder  bulk  gas  temperature 
'   oulk  gas    J 


V 


°D 


T  reference  temoerature  R 

o 

T,  boundary  layer  gas  temper-       R 
boundary        t  / 

layer 

T  -,-,  cylinder  wall  temperature       °R 

TSTR  non-dimensionalized  temper- 
ature (eqn  23) 

U  internal  energy  in-lb 


cylinder  volume  in^ 

3 


V  reference  cylinder  volume       in 

VBDC  cylinder  volume  at  piston       in^ 

bottom  dead  center  position 

VTDC  cylinder  volume  at  piston       in^ 

top  dead  center  position 

VSTR  non-dimensionalized  cylinder 

volume  (eqn  20) 

7/  total  cylinder  work  in-lb  f 

W„rtw,,,         work  done  during  compression    in-lb.. 
comP  stroke 

HL  work  done  during  expansion      in-lb  - 

exp  stroke  1 

ratio  of  specific  heats 
«VCv> 
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CHAPTER  I 
INTRODUCTION 

The  reciprocating  piston-cylinder  system  is  a  common 
and  relatively  simple  device  used  in  many  different  applica- 
tions in  today's  mechanical  world.   Each  day  some  type  of 
reciprocating  piston-cylinder  system  impacts  on  our  lives 
whether  it  be  through  a  reciprocating  automotive  engine  or 
through  some  other  type  of  reciprocating  pumping  device  that 
may  provide  another  personal  service.   An  additional  example 
of  this  type  of  device  is  a  reciprocating  cryogenic  expander, 
whose  purpose  is  to  produce  a  reduction  in  the  temperature 
of  the  gas  passing  through  it. 

In  all  of  these  types  of  systems  work  is  produced 
through  the  expansion  of  a  fluid  against  a  piston.   For  an 
ideal  expansion  engine  the  expansion  stage  should  be  isen- 
tropic  and  the  work  per  unit  mass  should  be  expressed  as: 

W  =  Cp  (T1  -  T2)  (1) 

However,  because  this  fluid  expansion  is  not  an  ideal  pro- 
cess, losses  and  inefficiencies  are  evident.   Some  examples 
of  typical  loss  mechanisms  usually  associated  with  through 
flow  expanders  are: 

1.   clearance  losses  -  Because  of  the  clearance 
volume  remaining  when  the  piston  is  at  top  dead  center,  the 
cylinder  cannot  fully  exhaust  the  entire  old  charge.   There- 


fore  the  fresh  charge  entering  the  piston  is  further  limited 
in  the  actual  volume  that  it  can  ideally  occupy.  The  clear- 
ance losses  then  define  these  variations  from  the  ideal 
case.   The  actual  factors  responsible  for  the  clearance  loss 
include  much  more  complex  processes  such  as  blow  in-blow  out 
loss  effects  and  the  intricate  influences  of  the  mass  flow 
rate.   These  processes  are  extremely  involved  and  their 
evaluation  here  is  not  required,  so  then  just  calling  atten- 
tion to  their  existence  will  suffice. 

2.   valve  losses  -  Because  the  inlet  and  exhaust 
valves  are  responsible  for  a  resistance  to  the  fluid  flow, 
there  will  be  a  pressure  drop  across  the  valves  that  will 
account  for  some  loss. 

3«   piston  ring  losses  -  There  will  be  some  loss 
associated  with  the  leakage  past  the  piston  rings.   Also 
any  friction  between  the  piston  rings  and  cylinder  walls 
will  generate  some  heat  that  will  be  transferred  to  the 
fluid  and  add  to  the  loss. 

The  above  examples  are  certainly  not  the  only  loss  mech- 
anisms evident  in  a  piston-cylinder  arrangement,  but  they 
are  usually  considered  the  prominent  factors. 

This  experimental  investigation  v/ill  consider  the  loss 
mechanisms  evident  in  a  closed  cylinder.   No  inlet  or  ex- 
haust processes  will  be  evaluated.   The  mass  of  the  working 
fluid  will  remain  constant  throughout  the  compression  expan- 
sion cycle  being  examined.   Therefore  the  causes  of  the 
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losses  evident  in  the  cylinder  will  simplify,  as  the  heat 
transfer  loss  will  be  the  only  prominent  loss  mechanism. 
This  will  then  permit  the  determination  of  the  actual  effect 
that  the  heat  transfer  loss  has  on  the  system. 

I. A.   Calculation  of  Loss 

To  evaluate  the  thermodynamic  processes  occurring  in 
the  closed  piston-cylinder  the  first  law  of  thermodynamics 
will  be  utilized.   In  its  most  basic  form  the  first  law 
states  that 

dG  5  dU  +  dl7.  (2) 

If  the  pressure  of  the  gas  is  spatially  uniform  at  each 
state  of  the  cycle  and  the  work  done  by  the  gas  is  re- 
stricted to  the  displacement  of  the  piston,  then  the  work 
term  can  be  defined  as 

dff  =  JpdV.  •  (3) 

For  an  ideal  gas  the  change  in  internal  energy  of  the  gas  is 

dU  =  mC   dT,  ,,  (L) 

v   -bulk  gas.  K^' 

The    first   lav/   can  then  be    expressed  as 

/Pdv  =  do  -  *cv  dTbulk  gas>       (5) 

If  the  process  is  assumed  to  be  ideal  and  adiabatic, 
then  the  heat  transfer  in  the  system  will  be  zero  (dO  =  0) 

1  1 


and  the  first  law  will  then  simplify  to 

/Pdv  =  mCv  dTbulk  gas>  (6) 

When  dT  is  eliminated  by  substitution  of  the  ideal  gas 
equation  of  state,  equation  (5)  integrates  to 

PV  =  constant.  (7) 

The  construction  of  the  trace  for  this  process  on  a  pres- 
sure vs.  volume  diagram  will  yield  a  plot  in  which  the  com- 
pression and  expansion  processes  share  the  same  line  and 
therefore  describe  no  enclosed  area  (Figure  1).   There  will 
be  zero  net  work  done  on  the  working  fluid. 

If  the  process  is  not  ideal,  the  compression-expansion 
trace  will  not  follow  the  PV  =  constant  path.   Plotting 
the  non-ideal  expansion  trace  on  a  pressure  vs.  volume 
diagram  will  yield  a  closed  plot  circumscribing  some  en- 
closed area  (Figure  2).   This  enclosed  area  represents  the 
net  work  done  on  the  working  fluid  during  the  cycle  and  also 
represents  the  total  losses  that  the  system  experiences 
during  the  cycle. 

Therefore  the  total  loss  per  cycle  in  the  piston- 
cylinder  arrangement  can  be  expressed  as  the  integral  of 
pressure  times  a  differential  volume  around  the  entire  cycle 

^PdV  =  loss.  (8) 

This  loss  is  specifically  associated  with  just  the  bulk  gas 
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in  the  cylinder.  The  total  loss  per  cycle  can  also  be  ex- 
pressed in  a  different  form  by  considering  the  loss  as  the 
integral  of  the  rate  of  heat  transfer  from  the  bulk  gas 
through  the  wall  around  the  entire  cycle 


Oq  dt  =  loss,  (9) 


1.3,   Boundary  Layer  Effects 

Adjacent  to  the  inside  of  the  cylinder  wall  a  small 
volume  of  gas  with  finite  mass  will  permanently  exist 
through  which  the  heat  flux  must  travel  in  order  to  reach 
the  cylinder  wall.   As  the  piston  moves  in  the  cylinder,  the 
depth  of  this  thin  layer  of  gas  will  vary  constantly  due  to 
well  defined  boundary  layer  effects.   If  the  cyclic  heat 
flow  entering  and  exiting  the  boundary  layer  is  considered, 
the  importance  of  the  boundary  layer  effect  can  be  illus- 
trated (Figure  3)» 

For  an  entire  cycle  the  total  loss  from  the  bulk  gas  to 
the  wall  (d)q  dt)  should  be  equal  to  the  loss  from  the  bulk 
gas  to  the  boundary  layer  (<fq-  dt)  and  should  also  be  equiv- 
alent to  the  loss  from  the  boundary  layer  to  the  wall 
((£q2  dt). 


(bq   dt  =  aq-  dt  =  mq-  dt  (10) 


Because  the  heat  capacitance  of  the  cylinder  wall  is  very 
great  when  compared  to  the  heat  capacitance  of  the  gas, 
the  wall  temperature  can  be  assumed  to  be  constant  when 

1Jf 
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BOUNDARY  LAYER 


TV  -,-,  =  CONSTANT 
'  wall 
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considering  a  single  cycle. 

If  the  small  volume  of  gas  comprising  the  boundary 
layer  had  no  heat  capacitance  then  at  any  instant  during  the 
cycle  the  rate  of  heat  transfer  from  the  bulk  gas  to  the 
boundary  layer  would  equal  the  rate  of  heat  transfer  from 
the  boundary  layer  to  the  cylinder  wall, 

(instantaneous,  zero  heat     \     f  i  1 ") 
capacitance  in  boundary  layer] 

However,  the  presence  of  this  boundary  layer  would  suggest 
that  it  would  hold  some  finite  value  of  heat  capacitance 
between  the  temperature  extremes  of  the  bulk  gas  and  the 
cylinder  wall.   Assuming  that  the  boundary  layer  does  sup- 
port some  heat  capacitance,  then  as  the  piston  moves  through 
the  cylinder  the  actual  rate  of  heat  transfer  from  the  bulk 
gas  to  the  boundary  layer  should  instantaneously  differ  from 
the  rate  of  heat  transfer  from  the  boundary  layer  to  the 
cylinder  wall. 

/    /      /instantaneous,  finite  heat   \     (\p\ 
•my   I  capacitance  in  boundary  layeri     ^   ' 

I.C.   Phase  Difference  in  Heat  Flux 

The  rate  of  heat  transfer  and  therefore  the  loss  expe- 
rienced from  the  system  depends  heavily  upon  the  temperature 
difference  between  the  working  gas  and  the  cylinder  wall  and 
upon  the  properties  of  the  working  gas. 

a-  =  hA  (Tgas  "  T.Jall)  (13) 
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During  both  the  compression  and  the  expansion  stroke  of  a 
cycle,  the  system  will  pass  through  a  point  where  the  gas 
temperature  will  equal  the  cylinder  wall  temperature.   At 
these  two  points  the  rate  of  heat  transfer  will  equal  zero 
as  the  temperature  difference  will  equal  zero  (T___  -  T,.rolT 

*  qcIo      •<  aj.  J. 

=  0),   Consider  a  system  in  which  the  cylinder  wall  boundary 
layer  has  zero  heat  capacitance.   Here  the  temperature  of 
the  gas  in  the  layer  adjacent  to  the  wall  will  be  approxi- 
mately equal  to  the  bulk  gas  temperature.   Therefore  at  the 
points  in  the  cycle  where  the  bulk  gas  temperature  equals 
the  cylinder  wall  temperature,  the  process  will  be  adiabatic 
and  the  rate  of  heat  transfer  will  equal  zero. 

However,  if  the  system  is  assumed  to  have  finite  heat 
capacitance  in  the  thermal  boundary  layer,  the  temperature 
of  the  gas  in  the  layer  adjacent  to  the  wall  could  differ 
significantly  from  the  bulk  gas  temperature.   In  this  system 
the  points  of  zero  rate  of  heat  transfer  will  not  occur 

where  T«  -,,  ^„_  =  Tmnll      but  will  occur  at  the  joints  where 
bulk  gas    wall, 

T  =  T 

boundary  layer    wall. 

To  illustrate  this  difference,  if  the  compression  ex- 
pansion process  is  plotted  on  a  pressure  vs.  volume  diagram, 
the  points  of  zero  rate  of  heat  transfer  will  occur  where 
the  process  is  adiabatic  at  that  point.  Previously  in  this 
section  it  was  shown  that  the  process  will  be  adiabatic  at 
the  point  -/here  the  trace  is  tangent  to  a  line  where  PV  = 
constant.   Figure  if  illustrates  a  pressure  vs.  volume 
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diagram  for  a  system  where  the  points  of  tangency  to  the 

adiabatic  lines  coincide  with  the  points  where  the  bulk  gas 

temperature  equals  the  cylinder  wall  temperature. 

If  a  system  with  finite  heat  capacitance  in  the 

boundary  layer  is  considered,  the  points  at  which  the  bulk 

gas  temperature  equal  the  cylinder  wall  temperature  will 

not  coincide  with  the  points  at  which  the  system  trace  is 

tangent  to  the  adiabatic  lines.   Figure  5  illustrates  this 

difference.   In  both  the  compression  and  expansion  stroke, 

a  phase  difference  is  identified  that  is  the  difference  (in 

time  or  degrees  of  rotation)  between  the  points  of  actual 

zero  rate  of  heat  transfer  in  the  system  (points  where  the 

process  is  perfectly  adiabatic  where  T,    j    -,  = 

J  boundary  layer 

T  ,,)  and  the  points  of  apparent  zero  rate  of  heat  transfer 
in  the  system  (points  where  Tbulk  0.as  =  T,vall), 

Therefore  the  effect  of  the  boundary  layer  on  the  heat 
transfer  in  the  cylinder  must  be  determined  in  order  to 
evaluate  the  impact  of  the  boundary  layer  on  the  total 
system  loss. 

Other  methods  are  certainly  available  that  would  facil- 
itate calculation  of  the  net  cyclic  loss  and  any  phase  dif- 
ference experienced  in  the  points  of  zero  rate  of  heat 
transfer.  One  of  the  most  direct  and  accurate  methods  re- 
quires the  plotting  of  a  temperature  vs.  entropy  diagram 
for  the  cycle.   Figure  6  is  an  example  of  such  a  trace.   The 
area  enclosed  by  the  curve  represents  the  net  work  done  on 
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the  working  fluid  which  is  also  equivalent  to  the  total  loss 
experienced  during  the  cycle.   For  any  single  cycle  the 
following  statement  should  be  true: 


c^PdV  =<^TdS  =  loss.  (IZf) 


On  this  temperature  vs.  entropy  diagram,  the  points  on  the 

trace  that  are  tangent  to  a  vertical  line  signify  points 

where  the  change  in  entropy  is  equal  to  zero.   At  these 

points  the  process  is  perfectly  adiabatic  and  signify  points 

of  zero  rate  of  heat  transfer.   The  temperature  of  the  gas 

adjacent  to  the  cylinder  wall  at  these  points  should  be 

equal  to  the  temperature  of  the  cylinder  walls.   If  these 

points  of  actual  zero  rate  of  heat  transfer  differ  from  the 

points  on  the  graph  of  apparent  zero  rate  of  heat  transfer 

(where  T,  -,,  _   =  T,„  , , ) ,  then  some  lohase  difference  exists, 
bulk  gas    wall  7 

In  Figure  6  a  phase  difference  is  evident  in  both  the  com- 
pression and  expansion  strokes. 

The  existence  of  this  phase  difference  must  be  proven 
and  then  evaluated  in  order  to  provide  information  necessary 
to  quantify  the  impact  of  the  thermal  boundary  layer  effects 
on  the  total  system  loss. 
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CHAPTER  II 
DESCRIPTION  OF  TEST  EQUIPMENT  AND  TESTING  PROCEDURE 

II. A.   Experimental  Expander 

The  apparatus  used  to  investigate  the  theoretical  pro- 
cesses of  the  expander  was  composed  primarily  of  the  motor 
and  cylinder  body  of  a  modified  air  compressor.   The  active 
cylinder  was  bolted  directly  on  top  of  the  original  com- 
pressor cylinder.   This  cylinder  included  the  gas  charging 
port,  temperature  and  pressure  measuring  instruments  and  had 
the  additional  capability  of  water  cooling  the  cylinder,   A 
diagram  of  this  apparatus  is  shown  in  Figure  7.   The  piston 
utilized  to  compress  and  expand  the  gas  was  housed  in  the 
active  cylinder.   The  active  cylinder  bore  was  2,00  inches 
and  the  piston  stroke  was  3»00  inches.   During  operation  the 
piston  remained  within  the  cylinder,  while  its  maximum  up- 
ward travel  at  the  top  dead  center  position  moved  the  piston 
to  the  upper  lip  of  the  cylinder. 

The  cylinder  extension  housed  a  stationary  piston  that 
formed  the  head  of  the  cylinder  volume.   A  threaded  bolt 
extending  through  the  top  of  the  extension  cylinder  and  into 
the  stationary  piston  provided  the  additional  capability  of 
changing  the  volume  ratio.   By  raising  or  lowering  the  sta- 
tionary piston  through  its  four  inch  range  with  this 
threaded  bolt,  the  cylinder  provided  volume  ratios  ranging 
from  1.5  to  2.4« 
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II. 3.   Pressure  Measurement 

A  strain  gage  pressure  transducer  ivas  utilized  to 
measure  cylinder  pressure.  The  pressure  measurement  is 
determined  in  this  type  of  transducer  by  the  fact  that  the 
resistance  of  the  transducer  strain  gages  on  the  pressure 
sensitive  surface  changes  in  proportion  to  the  pressure  on 
the  surface.   The  transducer  output  is  then  proportional  to 
the  cylinder  pressure.   The  transducer  mounted  in  the  cylin- 
der head  was  of  the  range  0  to  100  psi  and  was  calibrated 
at  if. 342  p si/millivolt  for  a  6.0  volt  input.   In  order  to 
calculate  the  cylinder  pressure,  the  transducer  output  at  a 
known  pressure  had  to  be  determined.   Therefore  prior  to 
recording  any  experimental  data,  the  transducer  input  at 
atmospheric  pressure  was  recorded  to  provide  the  required 
pressure  reference.   To  make  this  step  easier  a  balance  box 
was  installed  with  the  pressure  transducer  which  permitted 
the  setting  of  a  known  pressure  directly  on  the  oscilloscope. 

II. C.   Volume  Measurement 

The  method  utilized  to  measure  the  cylinder  volume 
consisted  of  a  circular  cam  mounted  eccentrically  on  the 
compressor  crankshaft.   Because  each  cycle  of  the  compressor 
piston  corresponded  to  one  revolution  of  the  compressor 
crankshaft,  the  angle  of  a  given  spot  on  the  circumference 
of  the  shaft  was  directly  proportional  to  the  position  of 
the  compressor  piston  and  therefore  was  also  proportional  to 
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the  cylinder  volume.   A  follower  was  mounted  on  the  outer 
diameter  of  the  circular  cam  so  that  the  displacement  of  the 
follower  v/as  proportional  to  the  displacement  of  the  piston. 
This  system  is  displayed  in  Figure  5. 

An  HP  7DCT  displacement  transducer  was  connected  to  the 
follower  and  provided  a  D.C.  output  that  was  proportional  to 
the  displacement  of  the  follower.   The  calibration  of  the 
transducer  provided  the  maximum  output  signal  when  the  cylin- 
der volume  was  maximum  at  the  piston  bottom  dead  center 
position,  3y  knowing  the  maximum  and  minimum  cylinder  volume 
for  the  experimental  run,  the  intermediate  values  for  cylin- 
der volume  could  easily  be  calculated  by  using  the  transducer 
output. 

II. D.   Cylinder  "Tall  Temperature  Measurement 

Four  thermocouples  were  utilized  to  measure  the  average 
cylinder  wall  temperature  during  testing.   Three  thermocou- 
ples monitored  the  temperature  at  different  locations  on  the 
outside  of  the  cylinder  wall.   The  remaining  thermocouple 
was  placed  on  the  outside  of  the  cylinder  cooling  water 
jacket.   The  output  of  the  thermocouples  was  displayed  using 
an  Autodata  Nine  Digital  data  logger  that  provided  a  contin- 
uous readout  in  degrees  Centigrade. 

II. 5.   Digital  Storage  Oscilloscope 

The  output  of  the  pressure  and  volume  transducers  v/as 
provided  as  input  to  a  Nicolet  Digital  Oscilloscope  Model 
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206.   The  signal  information  provided  to  the  scope  was  con- 
verted to  digital  form,  stored  in  a  buffer  memory  and  then 
was  displayed  graphically.   Digital  information  for  each 
point  remained  available  in  storage.   The  Nicolet  scope  has 
the  capability  of  storing  up  to  if096  data  points.   Because' 
both  pressure  and  volume  v/ere  being  monitored,  each  storage 
cycle  consisted  of  ZOUQ   pressure  points  and  20^-8  volume 
points.   The  scope  had  the  capability  of  graphic  display  of 
either  the  pressure  or  volume  trace  vs,  time  or  an  X-Y  axis 
plot  of  pressure  vs,  volume.   The  time  interval  between  data 
points  could  be  chosen  on  the  scope  to  ensure  that  an  entire 
cycle  of  data  would  be  available  in  a  single  storage  cell, 

II, F,   Testing  Procedure 

The  following  steps  were  performed  in  this  order  during 
the  data  acquisition  run: 

1,  The  position  of  the  stationary  piston  was  set 
properly, 

2,  The  piston  was  moved  to  bottom  dead  center 
position  and  then  checked  by  assuring  that 
the  output  of  the  volume  transducer  was  at  a 
maximum, 

3«  The  cylinder  was  opened  to  the  atmosphere  and 
the  pressure  transducer  output  was  set  on  the 
oscilloscope  and  recorded. 

A-,   The  cylinder  was  charged  with  gas  and  the 

pressure  transducer  output  recorded.   At  this 
point  must  ensure  that  the  cylinder  is  not 
charged  to  a  pressure  that  will  over-pressurize 
the  pressure  transducer  during  compression. 

5.   The  compressor  motor  was  turned  on  and  the 
oscilloscope  set  to  automatically  monitor  a 
single  cycle  of  the  cylinder, 
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6.  The  cylinder  wall  temperatures  were  recorded 

immediately  upon  completion  of  the  sample 
cycle. 

7.  The  compressor  motor  was  secured. 

8.  The  pressure  vs.  volume  data  was  monitored  on 
the  oscilloscope  to  ensure  that  the  trace  was 
satisfactory. 
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CHAPTER  III 
COMPUTER  REDUCTION  AND  INTERPRETATION  OF  DATA 

In  order  to  efficiently  interpret  and  display  the  raw 
data  extracted  from  the  Nicolet  oscilloscope,  the  Institute's 
CMS  computer  system  was  used.   The  data  was  entered  into  the 
computer  files  and  then  was  operated  on  to  provide  various 
thermodynamic  parameters  at  each  point.   Using  this  infor- 
mation thermodynamic  traces  for  the  cycle  were  plotted. 
Finally  the  losses  for  the  cycle  were  calculated  by  numer- 
ical integration  of  the  curves, 

III .A.   Computer  Inputs 

The  digital  pressure  vs,  volume  data  provided  by  the 
Nicolet  oscilloscope  was  input  into  a  data  file  that  was 
accessible  to  the  main  program.   Other  initial  inputs  re- 
quired included: 

1,  The  type  of  working  gas  used, 

2,  The  specific  gas  constant  for  the  gas  (R). 

3»   The  specific  heat  at  constant  pressure  for 
the  gas  (C  ). 

if.   The  specific  heat  at  constant  volume  for  the 
gas  (Cv). 

5.  The  cylinder  volume  at  the  top  dead  center 
position  (VTDC), 

6.  The  cylinder  volume  at  the  bottom  dead  center 
position  (VEDC), 

7.  The  total  number  of  data  points  (P,  V  pairs) 
entered  into  the  computer  data  file. 

8.  The  time  interval  between  data  ooints. 
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9.   The  piston  revolutions  per  minute. 

10,   The  temperature  readings  from  the  cylinder 
wall  thermocouples  in  degrees  Centigrade. 


III. 3.   Calculation  of  Reference  Points 

In  order  to  have  the  capability  of  accurately  comparing 
data  at  varying  volume  ratios,  compressor  speeds,  and  for 
different  gases,  computations  performed  on  the  data  were  done 
using  common  reference  points.   Therefore  a  standard  method 
was  required  to  determine  each  of  the  required  temperature, 
volume  and  pressure  reference  points. 

Cylinder  bulk  gas  temperature  was  the  first  reference 
point  determined  and  was  set  at  100  F.   This  temperature  was 
chosen  as  it  was  a  typical  average  cyclic  temperature  experi- 
enced in  previous  experimentation  done  with  the  same  experi- 
mental cylinder  apparatus.   When  expressed  on  the  absolute 
scale  the  reference  temperature  was  then  set  at  the  following 
value: 

TQ  =  559.67°  R 

The  method  used  to  determine  the  reference  volume  con- 
sisted of  calculating  the  geometric  mean  using  the  volume  of 
the  cylinder  at  the  top  dead  center  and  bottom  dead  center 
positions.   The  actual  formula  for  the  determination  of 
reference  volume  was 


V=(!f)T™c.  d5) 
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Because  the  cylinder  volume  was  expressed  in  units  of  cubic 
inches,  the  reference  volume  was  also  in  cubic  inches. 
The  determination  of  the  reference  pressure  point 
required  the  prior  determination  of  the  reference  volume 
point.   Using  the  reference  volume  point,  the  pressures 
corresponding  to  this  volume  for  both  the  compression  and 
expansion  strokes  were  found.   The  reference  pressure  was 
then  determined  to  be  the  average  of  the  compression  and 
expansion  pressure  values  at  the  reference  volume. 


P     +  p 
comT}    exiD 


P0  =  "uu"  g  ^u  (16) 


Reference  pressure  was  expressed  in  units  of  pounds  force 
per  square  inch. 

III.C.   System  Mass  Calculation 

In  order  to  calculate  the  bulk  gas  temperature  and 
entropy  at  each  data  point  the  system  mass  had  to  be  deter- 
mined.  This  calculation  was  performed  using  the  ideal  gas 
equation  of  state  and  t he  reference  temperature,  volume,  and 
pressure  values. 


P  V 

m  =  -^rr-2  (17) 

o 


The  system  mass  was  then  expressed  in  units  of  pounds  mass, 

III. P.   Calculation  of  Temperature  and  Entropy 

The  bulk  gas  temperature  at  each  data  point  was  then 

32 


determined  utilizing  the  ideal  gas  equation  of  state  and 
other  previously  determined  parameters. 


T.  = 
1 


R  m 


(18) 


The  temperature  values  were  expressed  on  the  absolute  scale 
in  degrees  Rankin e. 

The  entropy  value  for  each  data  point  was  calculated 
referenced  to  the  previously  determined  values  of  reference 
volume  and  pressure  by  the  relation 


S±   -  SQ   =  AS  =  raCv  In 


Pi 


"v ." 

+ 

mC^ 

In 

i 

V 

0 

(19) 


The  entropy  values  were  expressed  in  units  of  inch-pounds 
force  per  degree  Rankine. 

III. 5.   Calculation  of  System  Losses 

In  order  to  calculate  the  system  losses  during  the 
cycle,  cyclic  traces  of  pressure  vs.  volume  and  temperature 
vs.  entropy  were  plotted.   The  area  enclosed  within  either 
of  these  closed  curves  should  equal  the  total  system  loss. 


rfpdV  =  <^TdS  =  loss 


(1if) 


The  integration  was  performed  numerically  using  the  trap- 
ezoid rule  at  each  successive  data  point.   Total  area  under 
the  expansion  curve  was  considered  negative.   The  total 
system  loss  was  then  described  in  the  units  of  negative 
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inch-pounds  force. 

An  attempt  was  made  to  calculate  the  area  using  a  much 
more  involved  and  typically  more  accurate  five  point  Gaus- 
sian quadrature  integration  method.   However,  because  of 
the  scatter  of  the  digital  data  points  near  the  ends  of  the 
curves,  this  method  proved  no  more  accurate  than  the  much 
simpler  trapezoid  integration  method. 

III.F.   Non-Dimensional  Parameter  Calculations 

In  order  to  provide  the  additional  capability  of  com- 
paring the  computer  output  to  the  output  data  of  any  other 
similar  apparatus,  notwithstanding  differences  in  the  units 
involved,  non-dimensional  parameters  were  introduced.   The 
computer  routine  then  also  included  the  calculations  re- 
quired to  produce  the  output  of  non-dimensional  values  of 
pressure,  volume,  temperature  and  entropy.   Plots  were  also 
made  of  non-dimensional  pressure  vs.  volume  and  non-dimen- 
sional temperature  vs.  entropy  in  order  to  determine  the 
non-dimensional  loss  values. 

III.F. 1.   Non-Dimensional  Volume 

The  volume  data  points  7/ere  non-dimensionalized  by 
dividing  each  data  point  by  the  total  displaced  volume  of 
the  cylinder. 

VSTRi  =  (VBDC-VTDC)  '20) 
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III. P. 2.   lion-Dimensional  Pressure 

The  cylinder  pressure  data  points  were  no n- dim en si on- 
alized  by  dividing  each  data  point  by  a  mean  pressure  cal- 
culated from  the  compression  stroke.   This  mean  pressure  was 
defined  as  the  total  work  added  during  the  compression 
stroke  divided  by  the  total  displaced  volume  of  the  cylinder. 

VI 

P    -     "'comp  (7>    s 

mep  "  (VBDC-VTDC)  K       ' 

Each  pressure  data  point  was  then  divided  by  this  mean  pres- 
sure during  compression  to  find  the  non-dimensional  pressure 
values. 


p . 

PSTR.  =  p-±-  (22) 

meD 


III. P. 5.   Non-Dimensional  Temperature 

The  bulk  gas  temperature  data  points  were  non-dimension- 
alized  by  dividing  each  data  point  by  the  reference  temper- 
ature. 


Ti 
TSTP±  =  ^i  (23) 

o 


Ill.P.k.   Hon- Dimensional  Entropy 

The  entropy  values  were  non-dimensionalized  by  dividing 
the  calculated  entropy  values  by  the  mean  pressure  during 
compression  and  the  displaced  volume  of  the  cylinder  and 
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then  multiplying  by  the  reference  temperatur* 


Si  To 

S5TR  =  p    (VBDC-VTDC)  'l2if) 

mep 


III.F.5.   Non-Dimensional  System  Losses 

Non-dimensional  values  of  total  system  loss  were  cal- 
culated in  the  same  manner  as  the  dimensional  loss  values 
were  determined.   Flots  were  made  of  non-dimensional  pres- 
sure vs.  volume  and  non-dimensional  temperature  vs.  entropy 
and  the  enclosed  areas  were  calculated  by  numerical  inte- 
gration utilizing  the  trapezoid  method. 

III.G.   Computer  Output 

A  listing  of  the  Fortran  computer  program  is  displayed 
in  Appendix  A.  A  listing  of  the  test  data  file  used  in  the 
computations  is  presented  in  Appendix  B. 

The  output  of  the  computer  routine  includes  the  fol- 
lowing data: 

1.  Pressure,  volume,  bulk  gas  temperature  and 
entropy  at  each  data  point. 

2.  Non-dimensional  pressure,  volume,  bulk  gas 
temperature  and  entropy  at  each  data  point. 

3.  Piston  revolutions  per  minute. 

if.   Time  interval  between  data  points. 

5.  Reference  temperature  volume  and  pressure 
values. 

6.  Average  cylinder  wall  temperature. 

7.  System  mass. 
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8.  Pressure  vs.  volume  plot. 

9.  Temperature  vs.  volume  plot. 

10.  Temperature  vs.  entropy  plot. 

1 1 .  Non-dimensional  pressure  vs.  non-dimensional 
volume  plot. 

12.  Non-dimensional  temperature  vs.  non-dimen- 
sional volume  plot. 

13»   Non-dimensional  temperature  vs.  non-dimen- 
sional entropy  plot. 

llf.   Calculations  of  dimensional  and  non-dimen- 
sional system  loss. 

An  example  of  the  computer  output  is  displayed  in  Appendices 
C  and  D. 
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CHAPTER  IV 
RESULTS 

The  results  of  the  single  data  set  used  are  presented 
in  the  computer  output  in  Appendices  G  and  D.   Helium  was 
utilized  in  the  cylinder  as  t  he  working  gas  for  the  exper- 
iment. 

The  final  pressure  vs.  volume  and  non-dimensional 
pressure  vs.  non-dimensional  volume  plots  resulted  in 
extremely  smooth  traces.   Each  depicted  what  could  easily 
be  identified  as  a  conventional  pressure  vs.  volume  trace 
for  this  type  of  apparatus.   Basically  the  same  comment 
could  be  made  about  both  the  temperature  vs.  volume  and 
non-dimensional  temperature  vs.  non-dimensional  volume 
plots.   However,  a  greater  degree  of  scatter  seemed  to  be 
evident  in  the  temperature  vs.  volume  plots  than  in  the 
pressure  vs.  volume  plots. 

In  the  temperature  vs.  entropy  and  non-dimensional 
temperature  vs.  non-dimensional  entropy  plots  the  scatter 
is  much  more  evident.   The  shape  of  the  trace  also  is 
somewhat  different  from  a  shape  that  one  might  expect  the 
plot  to  have  after  first  seeing  the  pressure  vs.  volume  and 
temperature  vs.  volume  traces.   Figure  9  is  a  single  line 
curve  approximation  of  the  TSTR  vs.  SSTR  computer  trace 
from  Appendix  D.   The  figure  shows  the  adiabatic  and  iso- 
thermal limits  of  a  plot  of  this  type.   The  adiabatic  limit 
being  a  curve  where  the  compression  and  expansion  strokes 
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TSTR  VS.  SSTR  SMOOTH  CURVE  TRACE  WITH  ADIA3ATIC  AIID  ISOTHERMAL 
LIMITS  AND  COMPRESSION  AND  EXPANSION  STROKE  PHASE  DIFFERENCES 

FIGURE  9 
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approach  being  perfectly  adiabatic  while  the  isothermal 
limit  approaches  a  curve  where  the  compression  and  expansion 
strokes  are  perfectly  isothermal.   The  TST3  vs.  SSTR  plot 
has  been  drawn  on  a  graph  with  axes  of  equal  length  so  that 
the  curve  has  not  been  mistakenly  lengthened  in  either 
direction  to  give  the  false  appearance  of  being  more  adia- 
batic or  isothermal.   Nevertheless,  the  trace  seems  to 
approach  the  adiabatic  limit  in  the  high  temperature  portion 
of  the  plot,  while  approaching  the  isothermal  limit  in  the 
low  temperature  portion  of  the  plot.   On  the  compression 
stroke  the  trace  becomes  perfectly  adiabatic  at  a  point 
early  in  the  stroke  at  TST2  s  .365  (T  =  kSL°R) .      On  the 
expansion  stroke  the  trace  also  becomes  perfectly  adiabatic 
at  a  point  early  in  the  stroke  at  TSTS  »  1.15  (T  =  6H°R)  • 
At  each  of  these  two  points  that  actual  rate  of  heat 
transfer  in  the  system  through  the  cylinder  walls  equals 
zero. 

The  points  of  apparent  zero  rate  of  heat  transfer  occur 
at  T  =  T  -, 2  =  539. 2°3.   The  difference  between  these  points 
and  the  points  of  actual  zero  rate  of  heat  transfer  is  the 
phase  difference.   In  both  the  compression  and  expansion 
strokes  the  phase  difference  is  leading  as  the  point  of 
actual  zero  rate  of  heat  transfer  occurs  prior  to  reaching 
the  point  of  apparent  zero  rate  of  heat  transfer.   The 
actual  Dhase  differences  can  be  commuted  to  be: 


J+0 


compression  stroke   45    leading 
expansion  stroke     64    leading 


The  actual  system  loss  figures  computed  were: 

PdV  loss  =  -24.931  in-lb- 

TdS  loss  =  -25.107  in-lb^. 

I 

Non-dimensional  PdV  loss  =  -0.0584 
Non-dimensional  TdS  loss  =  -0,0588 

In  both  the  dimensional  and  non-dimensional  cases,  the  PdV 
and  TdS  losses  showed  less  than  one  oercent  difference. 
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CHAPTER  V 
DISCUSSION 

V.A.   Data  Point  Scatter 

The  scatter  evident  in  the  temperature  vs.  volume  plots 
and  particularly  in  the  temperature  vs.  entropy  plots 
deserves  explanation.   In  these  plots  this  scatter  occurs 
primarily  in  the  center  section  of  the  trace,  away  from  the 
areas  of  minimum  and  maximum  cylinder  volume.   At  the  upper 
and  lower  end  portions  of  all  the  plots,  the  data  points  are 
generally  packed  tightly  together  in  an  orderly  fashion. 
This  is  due  to  the  fact  that  at  the  ends  of  the  compression 
and  expansion  stroke,  the  piston  is  decelerating  rapidly, 
halting  abruptly  and  then  again  rapidly  accelerating  av/ay 
from  its  extreme  position.   This  stopping  and  starting 
motion  requires  a  great  deal  of  time  to  cover  a  small  dis- 
tance in  comparison  to  the  piston' s  much  more  rapid  pace 
at  the  mid-stroke  position  of  its  travel.   Because  the  data 
points  are  taken  at  equal  time  intervals,  many  more  data 
points  are  taken  at  the  extreme  positions  of  the  plots 
providing  the  tighter  packing  and  the  generally  smoother 
appearance  in  those  areas. 

Near  the  mid-stroke  positions  of  the  cycle  a  different 
situation  exists.   Here  the  piston  is  moving  much  more 
rapidly  relative  to  the  velocity  experienced  at  the  extreme 
cylinder  positions.   Again  as  the  data  points  are  taken  at 
equal  time  intervals,  a  greater  distance  should  exist 


between  points  on  each  trace  near  the  raid-stroke  position. 
This  situation  alone  does  not  account  for  the  scatter,  but 
along  with  two  other  possible  causes,  the  existence  of  the 
scatter  can  be  better  explained. 

First,  the  Nicolet  oscilloscope  samples  20 if 3  volume 
and  29A-8  pressure  data  points  per  scope  storage  cycle  at 
equal  time  intervals  between  each  point.   The  procedure  of 
the  scope  is  to  sample  one  volume  data  point,  wait  the 
appropriate  time  interval  and  then  sample  one  pressure  data 
point,  continuing  in  this  manner  until  the  storage  cycle  is 
full.   The  scope  does  not  have  the  capability  to  sample 
pressure  and  volume  simultaneously.   Therefore  the  resulting 
pressure-volume  data  point  pair  formed  by  matching  a  volume 
data  point  with  the  following  pressure  data  point  does  not 
exactly  coincide  with  that  pressure-volume  data  point  pair 
that  would  occur  if  the  capability  existed  to  sample  each 
pair  simultaneously.   This  fact  brings  a  small  amount  of 
error  into  the  position  of  each  data  point. 

Second,  the  digital  nature  of  the  data  generates  addi- 
tional error.   The  accuracy  of  the  digital  data  output  of 
the  Nicolet  was  limited  to  ,05  millivolt,   With  the  large 
number  of  data  points  taken,  the  changes  in  output  voltage 
between  many  points  was  less  than  .05  millivolt.   Therefore 
voltage  changes  of  less  than  .05  millivolt  were  rounded  off 
to  the  next  ,05  millivolt  increment  creating  some  error  in 
the  data.   This  might  have  not  been  critical  except  for  the 
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fact  that  the  volume  data  was  being  read  on  the  100  milli- 
volt scale  on  the  oscilloscope,  and  the  range  of  the  data 
covered  only  3«00  -  17.90  millivolts  (roughly  ]$%   of  the 
range  of  the  scale).   The  100  millivolt  scale  was  the  best 
scale  to  use  as  it  was  the  smallest  scale  available  on  the  ■ 
scope  that  would  accept  the  range  of  the  volume  data.   How- 
ever, a  ZfX  or  5X  signal  amplifier  should  have  been  connected 
in  the  volume  signal  wiring  to  utilize  the  full  range  of  the 
100  millivolt  scale.   This  would  have  alleviated  the  problem 
existing  due  to  the  .05  millivolt  accuracy  of  the  Micolet 
oscilloscope. 

•  These  limitations  of  the  digital  oscilloscope  in  con- 
junction with  the  magnitude  of  the  mid-stroke  piston  velocity 
producing  a  greater  distance  between  data  points  in  this 
area  of  the  plots,  are  surely  a  sufficient  basis  with  which 
to  attempt  to  justify  the  scatter. 

V.3.   Shape  of  Temperature  vs.  Entropy  Plots 

The  actual  shape  of  the  temperature  vs.  entropy  plots 
exhibits  a  few  variations  from  what  might  be  expected  after 
seeing  the  other  plots.   In  the  high  temperature  end  of  the 
plot,  the  curve  seems  to  approach  its  adiabatic  limit, 
especially  on  the  initial  portion  of  the  expansion  stroke 
(Figure  10).   The  explanation  for  this  is  fairly  simple. 
At  the  start  of  the  expansion  stroke,  the  temperature  of  the 
bulk  gas  and  the  temperature  of  the  gas  in  the  boundary 
layer  are  both  greater  than  the  cylinder  wall  temperature. 
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NON-DIMENSIONAL  TEMPERATURE  VS.  NON-DIMENSIONAL  ENTROPY  PLOT 
ANNOTATED  WITH  MINIMUM  AND  MAXIMUM  VALUES 


FIGURE  10 
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However,  the  temperature  of  the  bulk  gas  is  still  signifi- 
cantly greater  than  the  boundary  layer  gas  because  the  lower 
temperature  cylinder  wall  has  a  greater  effect  on  the  temper- 
ature of  the  gas  in  the  thin  boundary  layer.   As  the  expan- 
sion stroke  begins,  the  pressure  drops  and  the  bulk  gas 
temperature  drops  rapidly.   The  temperature  of  the  boundary 
layer  gas  drops  also,  but  it  drops  at  a  much  slower  rate  as 
it  started  the  expansion  stroke  at  a  significantly  lower 
temperature.   Because  the  temperature  difference  between  the 
boundary  layer  gas  and  the  cylinder  wall  is  not  great  for 
the  first  thirty  degrees  of  the  expansion  stroke,  very  little 
heat  is  transferred  out  of  the  system  and  the  process  ap- 
proaches its  adiabatic  limit. 

At  the  low  temperature  end  of  the  plot,  near  the  end 
of  the  expansion  stroke  the  curve  seems  to  approach  its 
isothermal  limit  (Figure  10),   The  explanation  for  this  is 
very  similar  to  that  for  the  near  adiabatic  case.   Near  the 
end  of  the  expansion  stroke  the  boundary  layer  gas  temper- 
ature has  fallen  somewhat  lower  than  the  cylinder  wall  tem- 
perature.  In  the  last  thirty  degrees  of  the  expansion  stroke 
the  volume  in  the  cylinder  increases,  but  it  increases  at  a 
relative  rate  much  less  than  the  rate  of  increase  near  the 
start  of  the  expansion  stroke.   Therefore  the  change  in  the 
bulk  gas  temperature  will  also  change  at  a  rate  much  less 
than  the  rate  of  change  at  the  start  of  the  expansion  stroke. 
Even  though  the  process  approaches  its  isothermal  limit, 
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there  is  still  a  sufficient  temperature  difference  between 
the  boundary  layer  gas  and  cylinder  wall  to  permit  signifi- 
cant heat  transfer  from  the  system. 

The  temperature  vs.  entropy  plot  also  has  what  seems  to 
be  an  inflection  near  the  mid-stroke  position  on  the  compres- 
sion stroke  (Figure  10).   At  this  inflection  there  is  an 
abrupt  change  in  the  slope  of  the  compression  trace.   Because 
of  its  proximity  to  the  mid-stroke  position  on  the  curve, 
this  inflection  may  be  due  to  the  change  of  piston  accelera- 
tion to  deceleration  at  mid-stroke.   No  such  glaring  inflec- 
tion is  evident  in  the  expansion  curve. 

7,C.   Identification  of  the  Adiabatic  Points  in  the  C"cle 

The  identification  of  the  two  points  in  the  cycle  that 
are  adiabatic  is  extremely  important  in  determining  the 
actual  phase  difference  between  the  heat  transfer  and  the 
bulk  gas  to  wall  temperature  difference.   It  is  expected  that 
both  the  compression  and  expansion  strokes  should  each  have 
one  adiabatic  point  where  Tboundapy  layep  =  Tv;all#   Ideally 
this  should  occur  at  the  points  of  maximum  and  minimum  values 
of  entropy.  However,  because  of  the  scatter  in  the  tempera- 
ture vs.  entropy  plots  it  was  difficult  to  envision  the 
actual  location  where  the  curve  became  tangent  to  the  ver- 
tical.  Therefore  to  identify  the  location  a  rough  curve  was 
passed  through  the  data  points  (Figure  9)  and  the  adiabatic 
points  were  then  evaluated  at  the  locations  where  the  curve 
passed  through  the  vertical. 
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To  check  these  adiabatic  locations,  lines  of  PV  = 

constant  (adiabatics)  were   constructed  on  the  pressure  vs. 
volume  plot  (Figure  11).   The  locations  on  the  trace  that 
were  tangent  to  an  adiabatic  line  would  signify  an  adiabatic 
point.   The  results  of  this  procedure  verified  the  selection 
of  the  adiabatic  points  in  Figure  9. 

V.D.   Effects  of  the  Boundary  Layer 

The  effects  of.  the  boundary  layer  on  the  rate  of  heat 
transfer,  or  more  specifically  on  the  losses  in  the  system 
are  most  important.   It  is  obvious  that  it  is  the  temperature 
of  the  boundary  layer  gas  and  not  the  bulk  gas  temperature 
which  is  the  driving  force  on  the  transfer  of  heat  through 
the  cylinder  wall.   If  the  boundary  layer  gas  temperature  is 
greater  than  the  temperature  of  the  cylinder  wall,  the  system 
will  lose  heat.   If  the  boundary  layer  gas  temperature  is 
less  than  the  temperature  of  the  cylinder  wall,  the  system 
will  gain  heat. 

Along  the  compression  stroke  and  part  of  the  expansion 
stroke  from  the  minimum  to  maximum  entropy  value  the  boundary 
layer  acts  to  cool  the  system  (Figure  10).   On  this  area  of 
the  curve  the  boundary  layer  gas  temperature  is  basically 
greater  than  the  cylinder  wall  temperature  and  heat  will  be 
transferred  out  of  the  system  through  the  cylinder  walls. 
The  heat  transfer  in  this  portion  of  the  curve  is  driven  by 
the  change  in  volume  and  not  the  change  in  the  bulk  gas  tem- 
perature of  the  system. 
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Along  the  expansion  stroke  from  the  minimum  entropy 
value  to  the  maximum  entropy  value  the  boundary  layer  acts 
to  heat  the  system  (Figure  10).   The  temperature  of  the  bulk 
gas  and  the  boundary  layer  is  falling;  however,  the  boundary 
layer  temperature  is  generally  less  than  the  cylinder  wall  . 
temperature  and  therefore  heat  will  be  transferred  into  the 
system  through  the  cylinder  walls. 

7.Z.   ZxTjlanation  of  the  Phase  Difference 

The  phase  difference  that  exists  in  the  zero  rate  of 
heat  transfer  location  is  a  very  important  effect  of  the 
boundary  layer  (Figure  9)»   The  explanation  of  the  cause  for 
the  phase  difference  should  be  fully  discussed. 

The  bulk  gas  temperature  is  an  important  parameter  in 
this  type  of  system.   Because  the  bulk  gas  temperature  is 
much  easier  to  measure  and  to  understand  than  the  more  dif- 
ficult to  measure  boundary  layer  temperature,  the  bulk  gas 
temperature  is  utilized  much  more  in  this  type  of  system. 
Therefore,  if  the  bulk  gas  temperature  is  utilized,  then 
when  deriving  the  apparent  expression  that  would  define  the 
rate  of  heat  transfer  from  the  system,  the  temperature  ex- 
tremes utilized  would  be  the  bulk  gas  temperature  and  the 
cylinder  wall  temperature  and  the  expression  would  be: 


a-  =  hA  (Tbulk  gas  -  -■/all' 


For  this  expression  the  points  of  zero  rate  of  heat  transfer 

are  the  joints  on  the  cycle  where  T^,,-,,      =  T   -, -, 

culk  gas    -wall. 
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The  existence  of  a  boundary  layer  with  finite  heat 
capacitance  alters  the  former  interpretation  of  the  problem. 
The  temperature  of  the  gas  in  the  layer  next  to  the  cylinder 
wall  will  not  equal  the  bulk  gas  temperature,  but  will  equal 
the  boundary  layer  gas  temperature.   Therefore,  the  expres- 
sion defining  the  rate  of  heat  transfer  to  the  cylinder  walls 
will  become: 

q  =  hA  (Tboundary  layer  "  TV;all)# 

This  expression  specifies  that  the  points  of  zero  rate  of 

heat  transfer  would  be  the  points  on  the  cycle  where 

T-u    j    i     =  T„  n    Because  the  bulk  gas  temperature 
boundary  layer    wall,  D 

is  not  equal  to  the  boundary  layer  gas  temperature  in  this 
system,  some  difference  will  exist  in  the  time  of  apparent 
zero  rate  of  heat  transfer  and  actual  zero  rate  of  heat 
transfer.   This  difference  is  the  zero  rate  of  heat  transfer 
phase  difference.   The  phase  difference  in  both  the  com- 
pression and  expansion  strokes  of  this  experimentation  was 
leading,  or  the  actual  zero  rate  of  heat  transfer  occurred 
in  the  cycle  before  the  apparent  zero  rate  of  heat  transfer 
(Figure  9). 

The  explanation  for  the  phase  difference  exists  in  the 
examination  of  the  temperature  history  of  the  boundary  layer. 
As  the  compression  stroke  is  started  the  boundary  layer  gas 
temperature  is  slightly  greater  than  the  bulk  gas  temperature. 
This  is  due  to  the  fact  that  the  cylinder  wall  temperature  is 
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greater  than  the  bulk  gas  temperature  and  the  cylinder  wall 
temperature  has  a  much  greater  effect  on  the  temperature  of 
the  boundary  layer  gas  rather  than  the  bulk  gas.   As  the  gas 
is  compressed  the  bulk  gas  and  boundary  layer  gas  temperatures 
rise,  but  because  the  boundary  layer  gas  temperature  started 
the  stroke  at  a  higher  temperature  than  the  bulk  gas,  the 
boundary  layer  gas  temperature  passes  through  the  cylinder 
wall  temperature  prior  to  the  bulk  gas  providing  a  leading 
phase  difference. 

As  the  expansion  stroke  is  started  both  the  bulk  gas 
and  boundary  layer  gas  are  cooled  by  the  drop  in  pressure. 
But  because  the  boundary  layer  gas  temperature  was  signifi- 
cantly lower  than  the  bulk  gas  temperature  at  the  beginning 
of  the  expansion  stroke,  the  boundary  layer  gas  temperature 
passes  through  the  cylinder  wall  temperature  prior  to  the 
bulk  gas  providing  a  leading  phase  difference. 
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CHAPTER  VI 
CONCLUSIONS 

One  of  the  original  objectives  of  the  work  included 
the  capability  of  direct  transmission  of  the  data  from  the 
experimental  expander  to  the  Institute's  CMS  computer  system. 
This  capability  required  the  use  of  a  converter  that  could 
properly  interface  with  the  Nicolet  storage  oscilloscope  and 
the  CMS  computer  system.   The  Electric  Power  Systems  Labora- 
tory was  assigned  the  task  of  designing  and  assembling  this 
interface  unit.   Due  to  manpower  shortages,  the  Electric 
Power  Systems  Laboratory  was  not  able  to  complete  the  as- 
sembly of  the  interface  unit  prior  to  the  submission  of  this 
paper.   The  unavailability  of  the  interface  unit  required  the 
data  to  be  read  point  by  point  digitally  from  the  Nicolet 
storage  oscilloscope  and  then  manually  entered  into  a  CMS 
computer  system  data  file. 

!.7ith  the  use  of  an  operable  interface  unit,  dozens  of 
cycle  runs  could  have  been  performed.   However,  the  lack  of 

the  interface  significantly  limited  the  number  of  cycle  runs 
that  could  be  evaluated.   Further  setbacks  occurred  when 
serious  problems  were  encountered  in  the  digital  and  graphic 
data  presentation  from  the  Nicolet  storage  oscilloscope.   The 
final  result  was  the  availability  of  only  one  good  set  of 
cycle  data  for  analysis. 

Although  the  technique  developed  in  this  paper  to 
analyze  the  data  was  used  on  only  one  set  of  experimental 
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data,  the  results  obtained  proved  that  the  analysis  method 
was  highly  credible,  7/hether  one  or  one  hundred  sets  of 
cycle  data  had  been  evaluated,  the  sane  method  would  have 
been  utilized. 

The  additional  runs  that  could  have  been  performed  with 
an  operable  interface  unit  would  have  greatly  increased  the 
scope  of  the  experiment.   Parameters  such  as  piston  RPM, 
Reynold's  number,  cylinder  volume  ratio,  cylinder  pressure 
and  type  of  working  gas  could  have  been  varied.   The  results 
that  would  have  been  presented  with  the  variations  in  the 
above  parameters  would  have  provided  insight  into  the  actual 
heat  transfer  losses  and  phase  difference  for  a  v/ide  range 
of  conditions.   This  general  information  could  then  be 
extrapolated  for  use  in  the  design  of  other  devices  of  this 
type. 
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APPENDIX  C 
CYCLE  PROGRAM  DATA  OUTPUT 
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